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The lossy propagation law (generalization of Lambert-Beer's law for classical radiation loss) for 
non-classical, dual-mode entangled states is derived from first principles, using an infinite-series of 
beam splitters to model continuous photon loss. This model is general enough to accommodate 
0^ \ stray-photon noise along the propagation, as well as amplitude attenuation. An explicit analytical 

. expression for the density matrix as a function of propagation distance is obtained for completely 

general input states with bounded photon number in each mode. The result is analyzed numerically 
for various examples of input states. For NOON state input, the loss of coherence and entanglement 
Q . is super exponential as predicted by a number of previous studies. However, for generic input states, 

1 where the coefficients are generated randomly, the decay of coherence is very different; in fact no 

' worse than the classical Beer-Lambert law. More surprisingly, there is a plateau at a mid-range 

interval in propagation distance where the loss is in fact sub-classical, following which it resumes 
the classical rate. The qualitative behavior of the decay of entanglement for two-mode propagation 
is also analyzed numerically for ensembles of random states using the behavior of negativity as a 
| function of propagation distance. 

O PACS numbers: 03.65.Ta, 06.20.Dk, 42.50.Lc, 42.50.St 
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We will derive the quantum analogue of the classical Beer-Lambert law for lossy propagation [l[ from first principles 
^ by modeling the propagation medium as a continuum series of linear optical scattering elements. The results we re- 
O" 1 port here suggest that the super-exponential propagation loss behavior of NOON states is highly special, and not likely 
to be shared by generic entangled states in the larger Hilbert space of the dual photon channel. The intuition that 
fSJ ■ entanglement embedded in a general state would decay similarly to the decay of entanglement in the lossy propagation 
of a NOON state appears to be faulty. Assuming even part of the coherent entanglement that survives during the 
propagation of a generic entangled state can be utilized to produce super-classical phase sensitivity using an appro- 
priate detection scheme Q , there appear to be many candidate states which are both robust against decoherence and 
non-classical enough to achieve significant advantage over classical light in an optimized quantum sensor architecture. 
This is consistent with and a generalization of previous results involving m-and-m states [H . 

OO ' 

The basic model is illustrated in Figure 1. Here we first analyze the propagation of single-mode quantum light. 
The model assumes a series of M identical beam splitters into which that the input modes are ao, d\, c?2, ■ ■ ■ , ^m, 
where ao is the incoming photon mode, and di , di, . . . , du are auxiliary channels possibly populated by stray light 
or thermal photons, but for the discussion here we will assume these channels have vacuum inputs. Similarly, the 
output modes are clm, s±, S2, • ■ ■ ,«Af, where a a/ is the output photon mode corresponding to the input channel ao, 
and Si , S2, ■ ■ ■ , sm are auxiliary channels modeling scattered and absorbed light along the propagation medium. The 
idea is to let the number of beam splitters M approach infinity in a controlled manner at the end of the calculation 
to extract the loss behavior of the quantum state propagating non-unitarily from the ao channel to the om channel. 

At each beam splitter k in Fig. 1, unitarity requires the relations 

flfe = r«n + Ld k , 

s k = La k -i + Td k , k = l, 2, ...,M, (1) 
where T and L are the complex transmission and reflection coefficients and satisfy the unitarity conditions 

|L| 2 + |T| 2 = 1 

LT + TL = (2) 
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FIG. 1: Diagram illustrating the series-of-beam-splitters model for propagation loss. 



In terms of the creation operators a\, d\, and s\, Eqs. (1) can be rewritten in the form 



4 = Ta l-\ 
4 = ^4-i 



■Ld\, 
Tdt, 



k = 1, 2, 



M 



(3) 



The input and output modes are connected by a unitary transformation that can be written in the form of a linear 
map 
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Using Eqs. (3), the (M + 1) x (M + 1) unitary matrix U can be written explicitly in terms of the reflection and 
transmission coefficients of the beam splitters: 
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It can be checked by straightforward calculation that the unitarity relations Eqs. (2) imply 
Since by Eqs. (4) — (5) 



aj = T M a\ { + Ls\ + LTs [ 2 + LT' z s\ 



(6) 
(7) 
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an incoming number Fock-state purely in the input mode 



(al) N \0) 



fNl 



is transformed to the outgoing state 
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fNl 



T M a} M + Ls\ + LTs\ + LT 2 s 
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10} 



(8) 



(9) 



As mentioned before, Eq. (8) assumes that the auxiliary input modes {di,d-2, - ■ ■ , d&r} are vacuum ports. More 
generally, our model can accommodate noise in the form of stray photons leaking into the propagation channel by 
simply replacing the input state Eq. (8) with a state of the form 



I^ n ) = -^(4) JV E c gi , 92 ,..., gM (4)«(4) 92 --(4) 9M |o> 

' 9l>92,- - ,<?Af 



To handle incoherent (such as thermal) input noise, one would have to enlarge the Hilbert space of the auxiliary input 
modes {d±, d%, ■ ■ ■ , g?a/ } to include their own environment modes, and trace over these secondary environment modes 
at the end of the calculation [4| . 

Going back to the input number state in the form Eq. (8), we can calculate the final output state (density matrix) 
in the out-mode (whose creation operator is a M ) by tracing over the loss (scattering) modes Si, S2, ■ ■ ■ , sm- 



Pout = Tr {si!S2i ._ SM} |Vw)(-0out| 
Substituting Eq. (9) for |Vw), this calculation gives 
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(10) 



(11) 



where, for ease of combinatorial manipulation, we renamed the output M-mode with index 0: so = o,m (equivalently, 
sj = a\#), and |fe) So = \k) SM . It is convenient to combinatorially manipulate Eq. (11) and rewrite it in the following 
form 



Pout 



E 



iV! 



N 



(JV-no)! 



, . v r:\no) s (n \s E i i 

no\{Jy - n )\ ^ ni\n 2 \...n M ' 

n =0 v ' n 1 ,...,n M =0 

J2f =1 n,=N -no 



(12) 



The continuum limit is defined by taking the number of beam splitters M to infinity, while keeping the "power 
lost" per unit length of propagation constant. Since each beam-splitter's contribution to power loss is given by \L\ 2 , 
mathematically this amounts to the limiting process (Fig. 2) 



Limit process for amplitude evolution: 
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Similarly, to preserve unitarity, a corresponding limit process must govern the evolution of the "phase" r\ of the 
complex transmission amplitude T: 



Limit process for phase evolution: 



(M — ► oo 
T= |T|e^ 
(b — > 



(14) 
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FIG. 2: The "Beer" 
propagation medium. 



limiting process where discrete beam splitters merge into a continuous 



Substituting Eqs. (13) and (14) into Eq. (12) yields our sought-for result for the x-dependent output density matrix 
resulting from an input Fock number state with N photons: 



N 



Pout (a:) = 



N 



c (l-e~H 



—fix \N —n 



n)(n\ . 



(15) 



The interpretation of Eq. (15) in terms of photon loss as a function of propagation distance is straightforward (see 
also [5] for an alternate derivation with a fixed amount of loss). 



Let's now turn to the analysis of dual-channel lossy propagation. We proceed similarly to the single-mode case, 
and adopt the notation for labeling the input, output, and scattering modes as illustrated by Fig. 3. Everything 
in the analysis above between Eqs. (1) and (7) proceeds independently for the two propagation channels a and 6, 
with corresponding series of beam splitters having in general distinct reflection and transmission characteristics. The 
dual-channel version of Eq. (7) is 
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As a first, important example, consider an input NOON state which has the form 



IVm) = j= (\N) a \0) b + \0) a \N) b ) = -i= ((4)^ + (bl) N ) |0) o |0) 6 . 
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According to Eq. (16), the output state is given by 

1 



|V>o 
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' T M h \ 
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As before, the output density matrix is obtained by tracing over the scattering modes: 



Pout 



Tr 



{si, ...,SM,tl, 



,t M }\lpout)(lpout\ 



(16) 



(17) 



(18) 



(19) 
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Dual channel entangled input 




FIG. 3: Modeling propagation loss for dual-channel entangled photon states with a series of 
beam splitters. 



Substituting Eq. (18) in Eq. (19) and carrying out some straightforward algebra, we reach the dual-channel analogue 
of Eq. (12), which we will suppress here to save space. Applying the Beer limit process (M — ► oo) defined by Eqs. (13) 
and (14) to this expression (whereby each channel a, b is associated with its own version of the extinction and phase- 
rotation coefficients fi a , fib, rj a , and rjb) yields the following continuum limit for the output density matrix of a NOON 
input state as a function of propagation distance x: 
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(20) 



Note that the coherence term is the entire expression on the second line of Eq. (20). This term is exponentially 
suppressed with an extinction coefficient given by iV/i, where fi = (fi a + /J,b)/2 is the average extinction coefficient 
of the two propagation channels. Hence the well-understood "super-exponential" loss of coherence (hence loss of 
entanglement) in the propagation of NOON states is once again verified [6|. 

After working through the NOON state example above, it is straightforward but rather cumbersome to extend our 
calculation to a completely general two-mode input state, the only restriction being that the maximum photon number 
TV in each mode a and b is finite. This most general input state can be written in the form 
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(21) 



Straightforward calculation in the spirit of the analysis presented thus far gives the result 
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where we also implemented a trivial generalization by allowing the output "point" to have different x coordinates x a 
and Xh along the two distinct propagation paths. Here and in what follows we adopt the combinatorial convention 
that the factorial of a negative integer is +00; thus the contributions to the above sum from, e.g., terms with I < p 
or p + I' < I vanish by virtue of the factorial terms in the denominator. A somewhat lengthy binomial-chase through 
the sums in Eq. (22) allows us to put it into the slightly more manageable alternative form: 

JV 

Pout{x a ,X b ) = E \p)a\q)b(p'\a(q'\b -X 

P; 1;V' ; ?'=0 

x e -^(P+p')f*aX a e ~i;(q+tl')LibXb e i(p~p')r)aX a e i(q-q')VbXb x 

X ai m ai+p' - P m +q' - q ( + p' - p)\ m\{m + q' - q)\\ 2 



X j~ ( l ~ P)K m ~ l) ] V p\q\p'\q'\ 

x (1 - e -^ x ") l -P (1 - e-" hXb ) m - q . (23) 



Before we begin analyzing the consequences of our main result Eq. (23) quantitatively, we mention one last general- 
ization which is quite straightforward. It is easy to incorporate x-dependent (variable) extinction and phase rotation 
coefficients p and rj into the above development merely via the substitutions in the limit process Eqs. (13)-(14) 



e - ft M (0 dc 

e iffvK)dC (24) 



With this generalization, the result Eq. (23) becomes 
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We now calculate the decay of coherence and entanglement numerically using the above formalism for NOON and 
generic (random) entangled states. The main results are illustrated by the plots in Figs. 4-7. 

All of our plots are generated assuming constant extinction and rotation coefficients fi — p a = Pb = 0.2km _1 and 
77 = 7] a = T]b = lkm - . The amount of coherence is calculated by computing the sum of absolute squares of the 
off-diagonal elements of the density matrix (the off-diagonal matrix norm squared) ; hence it can be interpreted as the 
coherence "power" that survives in the density matrix as it propagates through the lossy medium. Also note that in 
all three plots Figs. 4-7 the x axis is the propagation distance, and the y-axis depicts the logarithm of the dependent 
variable. 

Fig. 4 makes clear the "super-exponential" decay of coherence for a NOON state as a function of propagation distance. 
For the N = 10 NOON state, coherence power decays proportionally to exp(— 2Nfj,x) as expected from Eq. (20). 

In contrast to Fig. 4, Fig. 5 suggests an unexpected robustness of coherence in generic entangled states against 
propagation loss. Here we plot the coherence power as a function of propagation distance for two ensembles of 
entangled states, drawn randomly from two different probability distributions on the Hilbert space of all states of the 
form Eq. (21). One ensemble is drawn from a joint uniform distribution for the real and imaginary parts of ay, which 
is then normalized. The other ensemble is drawn from a true uniform distribution on the Bloch sphere in the complex 
space C N+1 . To generate the latter ensemble, we utilize the result that (zi, ■ ■ • z n )j + ■ ■ ■ + \z n \ 2 is a random 
vector on the unit sphere 5 12 ™" 1 in C™ with respect to the canonical volume form on S" 2 ™ -1 if Re(zj), lm(zj) are i.i.d. 
Gaussian random variables with zero mean 0,111 . We see from Fig. 5 that: 



• For a random entangled input state, coherence decays no faster than the classical loss rate. 



Log[coherence] 



-40 



■x(km) 



FIG. 4: Decay behavior of coherence power (the sum of absolute squares of the off-diagonal 
elements) for a N = 10 NOON state. 
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FIG. 5: Decay of coherence power for two ensembles of random entangled states (Eq. (21)) 
with N = 10. Note the slight plateau in decay at medium propagation distances. The plateau 
reverts back to exponential decay at large propagation distances; however, the overall decay 
remains no worse than that of classical light. The upper ensemble is drawn from a joint uni- 
form distribution for the real and imaginary parts of c%-, which is then normalized. The lower 
ensemble is drawn from a true uniform distribution on the Bloch sphere in the complex space 



where here N = 10. 
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• There is a mid-range plateau in loss, where decay is sub-classical, after which classical decay-rate resumes. This 
feature is present for all generic entangled inputs. 

• How much of the surviving generic coherence is useful for super-classical resolution will depend more on the 
survival rate of entanglement than on that of coherence. 



Log[negatiyity] 




FIG. 6: Behavior of negativity for the N = 10 NOON state and a random ensemble of entan- 
gled states (iV = 10). The ensemble is drawn from a true uniform distribution on the Bloch 
sphere in the complex space C JV+1 as in Fig. 5. 

To investigate the decay of entanglement in the density matrix that results from the propagation of a generic input 
state of the form Eq. (21), we use the concept of "negativity" [1, Consider a mixed state of a general bi-partite 
system with Hilbert space "Ka <8) 3~Cb given by 

P = E pij,kl\tAjB) &> (WbI , (26) 

i,j,k,l 

where {\iAjs}} is a separable orthonormal basis for %a ® 'Kb- The "partial transpose" of p with respect to the 
subsystem A is defined as the operator 

p tA = Yl Pkj,u\iAjB) ® (k A l B \ ■ (27) 

i,j,k,l 

Although the partial transpose p tA is symmetric and has unit trace, it is in general not a density matrix since it is 
not necessarily a positive operator. In fact, the separability of p is a sufficient (though not necessary) condition for 
the positivity of p tA . In other words, the negativity of p tA guarantees entanglement, but it is not necessary that p tA 
is negative for entanglement to be present. This makes negativity a partial entanglement measure which is relatively 
easy to compute. Quantitatively, negativity is defined to be the (non-negative) quantity 

E i A * A i' ( 28 ) 

i:\\ A < 

where the sum is over all eigenvalues of p tA which are negative. 
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In Fig. 6 we illustrate the behavior of negativity as a function of propagation distance for a NOON state and for a 
random ensemble of entangled states drawn randomly from the uniform probability distribution on the Bloch sphere 
in the Hilbert space of all states of the form Eq. (21). (See [ll[ for a different look at the behavior of entanglement 
under loss.) 

It is apparent that the "super-exponential" propagation loss behavior of NOON states is a highly special property not 
shared by more general two- mode entangled non-classical states of light. Whether the latent robust entanglement of 
generic two-mode states illustrated in Fig. 6 can be used to reach super-classical resolution in metrology applications 
depends on what detection schemes can be deployed on the output density matrix |l2j . as well as, from a practical 
point of view, on whether the generic states can be created reproducibly. These are questions we will investigate in a 
forthcoming paper [4(. 
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